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Abstract. We describe the spectral properties of underdoped cuprates as resulting from a momentum- 
dependent pseudogap in the normal state spectrum. Such a model accounts, within a BCS approach, for 
the doping dependence of the critical temperature and for the two-parameter leading-edge shift observed 
in the cuprates. By introducing a phenomenological temperature dependence of the pseudogap, which finds 
a natural interpretation within the stripe quantum-critical-point scenario for high-T c superconductors, we 
reproduce also the T c — T* bifurcation near optimum doping. Finally, we briefly discuss the different role 
of the gap and the pseudogap in determining the spectral and thermodynamical properties of the model 
at low temperatures. 

PACS. 74.25.Dw Superconductivity phase diagrams - 71.10.Hf Non-Fermi-Liquid ground states, electron 
phase diagram and phase transitions in model systems - 74.20.Fg BCS theory and its developments 

1 Introduction the critical temperature T c up to a crossover tempera- 
ture T* [Q. While T c increases with increasing doping, 

The non-Fermi-liquid behaviour of the normal phase of T * starts from much hi g ner values and decreases. The 

the cuprates has two major features: (i) nearby optimum tw0 temperatures merge around or slightly above opti- 

doping the in-plane resistivity is linear in T, signalling the mum doping. Angle-resolved photoemission spectroscopy 

absence of any other energy scale besides the temperature; (ARPES) experiments indicate that the pseudogap and 

(ii) in the underdoped regime photoemission and tunnel- the superconducting gap have the same momentum depen- 

ing experiments show that a pseudogap persists well above dence across T c , which almost resembles a ^2^2 -wave, 
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namely Z\k = Ao(<j)) cos(2</>), with = arctan(fc a / k x ) , is made however more involved by the presence of the 
and that both are tied to the underlying Fermi surface Landau renormalization factors [Q . 



In the BCS cZ-wave approach Aq is <p- independent 
and is proportional to T c . Here instead Ao((f>) is angle de- 
pendent and the ARPES spectra of the pseudogap state 
evolve smoothly and continuously in the superconduct- 
ing ones across the critical temperature. A smooth evo- 
lution is also observed in tunneling spectra ||. More- 
over, the pseudogap at different k points opens at dif- 
ferent temperatures. At T* a leading-edge shift {LE) ap- 
pears in the ARPES spectra of underdoped Bi2212 for 
momenta near the M = (±7r, 0); (0, ±7r) points of the Bril- 
louin zone. The LE is observed as a finite minimum dis- 
tance of the quasiparticle peak from the Fermi level in the 
non superconducting state, and leaves disconnected arcs 
of Fermi surface. When the temperature is lowered the 
LE regions around the M points enlarge and the arcs of 
Fermi surface reduce and shrink towards the nodal points 
of the corresponding d-wave superconducting gap below 
T c . At the same time the doping dependence of the mo- 
mentum structure of the LE is not trivial. As the doping 
is increased, the zero temperature LE at the M points, 
z4o(0), remains constant or decreases while the 

LE around the nodal points Aq(tt/A) seems to increase 



Many theoretical models have been proposed to obtain 
a non- Fermi-liquid behaviour and to describe a pseudogap 
state. A firm result is however that above one dimension 
the Landau Fermi-liquid theory is generically stable and a 
strongly singular effective potential is required to disrupt 
it 0. This result, together with the above phenomenol- 
ogy, suggests that a consistent description of the cuprates 
requires a strong momentum-, doping-, and temperature- 
dependent effective interaction. This interaction should af- 
fect the states near the M points of the Brillouin zone 
more effectively than along the diagonals and should have 
the temperature as the only energy scale around optimum 
doping. It was shown that in strongly correlated systems 
in the presence of additional attractive interactions (e.g. 
the Hubbard-Holstein model) and of long-range Coulomb 
forces, the exchange of quasi-critical charge (and spin) 
fluctuations provides such an effective electron-electron in- 
teraction both in the particle-particle and in the particle- 
hole channel |)|. Non- Fermi- liquid behaviour and strong 
pairing mechanism have in this way a common origin. 
These fluctuations arise near a finite-temperature instabil- 



ity line Tcdw(S) for charge-density wave or stripe-phase 

plffll and follow the rising of the critical temperature. r . , . , . . .... ,r^^r.\ 

L - FLJ formation, which ends m a quantum critical point (CJCP) 



Penetration depth measurements of the superfluid density 
p s (T) at low temperature probe the low-energy excitations 
around the nodal points in a d-wave superconductor, and 
therefore Aq{-k/4). The correspondence between ARPES 



at T = and S = S c near optimum doping P,^0|. As 
shown in Ref. ||, the effective electron-electron interac- 
tion near the charge instability has the form 



measurements of A (ir/4) and the slope of p s (T) at T = T . , . - V , . 

l/ e //(q,u;J ~ U 5— | 12 — : — , (1) 

«r + q — q c 2 - 17a; 



L. Benfatto, S. Caprara, and C. Di Castro: Gap and pseudogap evolution ... 3 

both in the particle-hole and particle-particle channels, particle-particle or particle-hole channel) are simpler to 

Here q = (q x ,q y ) and u> are the exchanged momenta and analyze and each of them shows relevant aspects of the 

frequencies in the quasiparticle scattering, U is a resid- physics of the cuprates. The interplay of the two channels 

ual repulsion, V is the strength of the attractive effec- is an open problem still under investigation. 

five potential, q c is the critical wave- vector related to the » c , . , , , , 

r ' ^ A first possibility is that the pseudogap opens due to 

charge ordering periodicity (q c — 2ir/\ c ). For physically . , , . . ,. , , iij 

° ° * ' ' v J J incoherent paring m the particle-particle channel, lead- 

relevant values of the parameters of the Hubbard-Holstein . , ^ . , , , , , . 

mg to a state where Cooper pairs around the M points 

model q c turns out to be (±0.28, ±0.86) or equivalently r j , rp. v, ^ , n n , , 

^ are formed at T > ±cdw with strong long- wavelength 

(±0.86, ±0.28) 101. In this case, therefore, q c connects the a , ,. D , , , , , . , 

v ' ' u ' ' ^ fluctuations. Phase coherence, which characterizes a real 

two branches of the FS around the M points and strongly n . , , , , 

1 ° J superconducting state, is established at a lower tempera- 

affects these states. The mass term k 2 = £~ 2 is the inverse mi ,i j_-rr ri , . . 

^ c ture i c , by coupling to the stiffness of the pairing near the 

square of the correlation length of the charge order and noc [ a i points |jTl| 
provides a measure of the distance from criticality. This 

In this paper we elaborate the other possibility, that 

is given by 8 — S c in the overdoped region, by T in the 

the transition to the superconducting state takes place in 

quantum critical region around S c and by T — Tcdw{5) 

the presence of a normal-state pseudogap parameter A p 

in the underdoped region, where Tcdw($) sets in a new 

resulting from interactions in the particle-hole channel. 

doping-dependent energy scale closely followed by T*(S). 

The issue arises of the interplay between the preformed 

The characteristic time scale of the critical fluctuations 

pseudogap in the p-h channel and the additional pairing 

is 7. The presence of a weak momentum-independent re- 
in the p-p channel. Having included most of the anomalous 

pulsion U together with a strong attraction of the order 

effects in the pseudogap formation, we determine T c via 

of —V/k 2 in the particle-particle channel (cfr. Eq. (|f|) fa- 

the BCS approach for the pairing in the p-p channel. Our 

vors d-wave superconductivity approaching optimum dop- 

model originates as a simple schematization of a system 

ing from the overdoped regime, within direct BCS calcu- 

interacting via the singular effective interaction (|]J) and 

lations JTo| . In the underdoped regime we expect that pre- 

is inspired to a similar model proposed by Nozieres and 

cursor effects of charge ordering are relevant to the pseu- 

Pistolcsi [l^] , with the inclusion of some specific aspects of 

dogap formation and extend up to a temperature Tq(S), 

the phenomenology of the cuprates. In Section 2 we discuss 

(the mean field temperature for CDW formation) higher 

the model for the normal-state spectrum in presence of a 

then T C dw ~ T* . 

pseudogap which has a d-wave form with amplitude A p . 
The two limiting cases when these precursors dominate Assuming at the beginning a costant A p , we discuss in 
the pseudogap formation in a single channel only (either Section 3 the general properties of our model, devoting 
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a particular attention to the doping and/or temperature the normal-state spectrum as 
dependence of T c , of the LE and of the superfluid density. c , I 2 , ai 2 

In Section 4 we introduce a modulation for A p , to take 

where 77 = +1(— 1) in the conduction (valence) band, et = 

care of the (5-dependence of the new energy scale set by 

— 2t(cosfc x + cosfcy) is the tight-binding dispersion law in 

the Tcdw(5) in the underdoped case. We assume that the 

the conventional metallic state (i.e. at A p — 0), A p is the 

pseudogap opens nearby a mean-field temperature Tq(S) 

pseudogap parameter, and /i is the chemical potential. A v 

for the onset of CDW. Tq(S) should follow the doping 

governs together with the chemical potential the LE open- 
dependence of Tcdw{$) in the underdoped regime and 

ing at the M points. The modulation of the pseudogap is 

produce a variation in the density of states, as revealed by 

given by the absolute value of 7k = (cosfcy — cosk y )/2, 

NMR and resistivity measurements on several compounds 

which vanishes along the diagonals of the Brillouin zone. 

[H. By a suitable fitting of Tq(S), we give at the end a 

Since ek vanishes along the lines k y — ±(n ± k x ), the 

phcnomenological description of the phase diagram of the 

two bands merge at the points kp = (±7r/2, ±7r/2). In 

cuprates, together with some physical quantities like the 

the undoped system the LE has the k-modulation given 

superfluid density, the specific heat and the leading edge. 

by the factor 7k. As the system is doped with respect to 
half-filling, a Fermi surface appears in the form of small 
2 The model pockets around the kp points. The LE, which vanishes 

along the Fermi surface, evolves thus continuously with 
Within the above scenario, we describe the pseudogap in respect to the undoped case. This is in contrast with the 
the normal state by means of a simplified model where a case of a k-independent pseudogap A p , which leads, in the 
k-dependent separation is present between a valence band lattice model, to a large Fermi surface in the doped sys- 
and a conduction band, as a result of a k-dependent ef- tern, with the chemical potential and the LE which evolve 
fective interaction in the particle-hole channel. Differently discontinuously with respect to half-filling, 
from Ref. |Q, we adopt a lattice electron model and as- We are aware of the fact that the complicated band 
sume that the pseudogap vanishes at some points of the structure of the cuprates is not simply fitted by the form 
Brillouin zone. Being interested on very qualitative as- (|J). For instance, the band structure can be improved 
pects of the evolution of the pseudogap state, we shall by extending the tight-binding expression beyond nearest 
mainly concentrate on a two-dimensional system related neighbours. However, this does not represent a severe lim- 
to the Cu02 planes, the third dimension being relevant itation to our approach, the generalization being straight - 
to establish the nature of the true transition and to cut forward. Moreover, the presence of the pseudogap leads, 
off the corresponding fluctuations. Accordingly, we model within our model, to the appearance of small pockets in 
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Fig. 1. Density of states at half-filling in the presence and in 
the absence of pseudogap for a hopping parameter t = 0.25 
eV. Observe in the case A p — 160 meV the splitting of the 
van Hove singularity which gives rise to two peaks at energies 
£ w ±A P , with A p < 21 

the weakly doped system. This issue is experimentally con- 
troversial, and it seems definitely confirmed that the two 
branches of each pocket, even when observed are not 
equivalent Jl4| . However, the presence of pockets around 
the kp points, which result from the oversimplification of 
our description, is not essential to reproduce the finite LE 
observed in the single-particle ARPES spectra near the M 
points, and does not play any significant role in the forth- 
coming analysis. Therefore, we shall not further discuss 
this aspect. 

The density of states corresponding to the band struc- 
ture (||) vanishes only at £ = — /i, and is finite elsewhere, 
as shown in Fig. 1. The van Hove singularity which exists 
at £ = — \x for A p = is thus split into two singularities 
separated by 4tA p /yJ(A p /2) 2 + 4t 2 . For A p < 2t, as we 
shall assume in most of our calculation, the energy range 
where the density of states in suppressed is of order 2A P . 



Having included most of the effects of the anisotropic 
potential in the k-dependence of the pseudogap Z\ p |7k|, 
for simplicity we assume that the onset of superconduc- 
tivity is produced within a BCS approach by a constant 
pairing interaction among the carriers in the Cooper chan- 
nel, apart again from a d-wave modulation. The relevance 
of the superconductive fluctuations of the phase of the or- 
der parameter will be discussed below, in parallel with the 
analysis of the properties of the superfluid density. 

We introduce therefore an intraband interaction term 
between time-reversed states 

d 2 k 



Hi 



d 2 k' 



(27T) 



2 C77,— k'XC77,k't j 



(3) 



where V is the pairing strength, Q is the volume of the 
system, c^" k g is the creation operator of the electrons in 
the ?7-band, and the factors 7k make the <i-wave sym- 
metry explicit. Introducing the order parameters A n = 
^(7kC^~ k jC^_ k |) we obtain the self-consistency equations 



A - - 
" ~ 2 ./ (2ir) 2 E nk 



^ * tanh 



1 



It, 



d 2 k 
(27P 



V 2 



A v (r, = ±1X4) 
, (5) 



where E,^ = y C r ^k + is the quasiparticle energy 

in the superconducting state. The last equation fixes the 
chemical potential \i for any given doping S with respect 
to half-filling. As most of the cuprates become supercon- 
ducting by doping with holes, we study the hole-doped 
regime S > 0, where the chemical potential at T = 0, in 
the absence of pairing, falls in the valence band. There- 
fore, at weak-coupling, we find the solution A^-i = A, 



6 L. Benfatto, S. Caprara, and C. Di Castro: Gap and pseudogap evolution ... 

A v —±i = for Eq. (^). We point out that, due to the conductivity in the undoped system, due to the fact that 
assumed band structure (Q), the spectrum is particle- hole at half- filling the density of states vanishes at the Fermi 
symmetric, and doping with electrons leads to the same energy. In our model, due to the d-wave nature of the pre- 
self-consistent solution, provided the role of the two bands formed gap, the value of V c is less than that found in Ref. 
is interchanged. [[l2) in presence of a constant gap. The equation that de- 

Adding to the Hamiltonian (|^) an interband interac- fines V c at T = 0, S = as a function of the pseudogap 
tion term, as discussed in Ref. fl2|| , induces the A v =+i amplitude A p is 

to be different from zero even at weak-coupling. The two 2t _ 1 f d 2 k 7^ 

Vc ~ 2 J (2ir) 2 ^(cos k x + cos kyY + (A p /2t)*T* 

order parameters are generically different, contrary to the 

so that V c « n 2 A p /2, for A p » 2t, and V c ~ 2t/\n(4t/A p ) 

assumption of Ref. |jl2|| . Indeed, the solution with one or- 

for Ap <C 2t. The former limiting case is however not re- 

der parameter (A^+i = A rt ——x) has a higher free energy, 

alistic: henceforth we shall definitely assume A p <C 2t. 

which may become equal under very specific conditions. 

Moreover, in the present scenario, the pseudogap is asso- 
ciated with scattering in the particle-hole channel. Pairing 

3 General properties of the model 

must therefore be absent at half-filling and we assume as 
In this section we shall analyze the main outcomes of our the relevant regime for the cuprates V < V c . 
model and their dependence on doping i5 and on the tern- Critical temperature Within the BCS approach, the 
perature T. The parameters A p and V at the beginning variation of the critical temperature with doping is con- 
will be assumed to be temperature and doping indepen- trolled by the density of states at the Fermi level. At 
dent, and not related to each other. This preliminary study A p — 0, we recover the spectrum of a normal metal with 
allows us to point out the consequences of the assumed k- nearest-neighbours hopping: at half-filling the Fermi en- 
dependent normal-state pseudogap, and to single out the ergy is at the van Hove singularity, and T c is maximum, 
regime of parameters suitable for a description of the phe- By increasing doping, the density of states at the Fermi en- 
nomenology of the cuprates. ergy decreases, and so does T c . This behaviour is excluded 

By solving the Eqs. (^)-(|^) we obtain the temperature by the experiments. On the contrary, when A p > 0, the 
and doping dependence of the superconducting gap A and critical temperature is zero at half-filling (if V < V c ). At 
of the chemical potential fi, and we determine the phase fixed finite doping, the critical temperature increases with 
diagram of the model in different physical regimes. increasing the pseudogap parameter A p , reaches a max- 

Before analyzing the two regimes of weak or strong imum, and then decreases. By doping, the critical tem- 
coupling at various doping, we investigate the existence perature follows the evolution of the density of states at 
of a finite strength V c of the interaction to have super- the Fermi energy, reaching a maximum when the Fermi 
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energy passes through the boundary of the pseudogap re- ARPES, penetration-depth and Andreev experiments 
gion in the density of states, i.e. at the doping S op t such The amplitude at the M points (LEm) is given 

that fi ~ —A p . For <5 > 8 op t, T c decreases with increasing essentially by the value of the LE in the normal state, i.e. 



doping. The resulting bell-shaped T c vs S curve captures 
the main features of the experimental results, as it will 
be more accurately presented in Section (4) with a doping 
dependent A p . 

The leading edge The quasiparticle spectrum, both in 
the normal and in the superconducting state, is charac- 
terized by a leading-edge shift LE, which is defined as 
LE((f) — nihi^fe .E^k, where k = /c(cos <p, sin0) []. Due 
to the point symmetries of the system our analysis can 
be limited to the quadrant < <f> < 45°. In the strong 
coupling regime the superconducting gap A 3> A p , so 
the doping variation of A dominates the evolution of the 
zero-temperature LE at the M points, which increases ap- 
proaching optimum doping, following T c , in contrast with 
the behaviour observed in the experiments in the under- 
doped regime. In the weak-coupling regime instead the 
superconducting gap is much smaller than the pseudo- 
gap. For S < S op t the LE of the superconducting state is 
controlled by two independent parameters (see Fig. (2)), 
whose existence is also suggested by the combination of 



LEm = 



2tA r 



V(4>/2) 2 +4* 2 



\n\ ^A p - \y\ (A p < 2t) 

(7) 



independent of T c , whereas the slope at the nodes is va = 
Aa((p — 7r/4) cx A cx T c , in agreement with the experi- 
mental finding [j§|,[|. Note instead that LE M ~ A for 
Ap = 0. 




10 15 20 25 30 35 40 45 
<p (degrees) 

Fig. 2. Zero-temperature LE vs <j> for t = 250 meV, 8 — 0.1, 
in the weak-coupling regime V = 87 meV, in the presence 
(x) and in the absence (+) of the normal-state pseudogap A v . 
Observe, for A p = 160 meV, the crossover around <f> ~ 12° to 
a regime dominated by the normal state pseudogap. 



The variation of the chemical potential accounts now 
for the main dependence of the LE on doping. By in- 



1 Within our mean-field treatment the excitations have a 
well defined energy and this quantity should be more cor- 
rectly called "one-particle excitation gap". However, by using creasing doping the normal-state LE M , controlled by the 
the term LE we want to put more emphasis on the distribu- pseudogap parameter A p according to Eq. @), decreases 
tion of spectral weight than on its coherent character, which and vanishes at a doping S c . If A p is doping independent 
is an artifact of our description of the cuprates in terms of a (as we are assuming in this preliminary discussion), the 



semiconducting band structure. 



S c for which LEm = coincides with the doping S opt for 
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the highest T c . At S > 5 C the LEm in the superconduct- 
ing state is then given only by the superconducting order 
parameter A, as shown in Fig. 3. At fixed doping, the 
temperature variation of the LE depends on the closing 
of the superconducting gap at T c , and eventually on the 
temperature dependence of A p , if any (see Fig. 4). 



0.035 




10 15 20 25 30 35 40 45 
<p (degrees) 

Fig. 3. Doping dependence of the zero temperature LE for 
t = 250 meV, A p = 160 meV and V = 87 meV. The 
LEm = LE(4> = 0) decreases with increasing doping, due to 
the reduction of the normal-state LEm given by Eq. (g). At 
S > S opt ~ 0.2 the normal state LE closes and LE(4>, T = 0) 
is that of a standard d-wave superconductor (see also Fig. 2). 
Instead, the slope of the LE at the node va oc A increases be- 
fore reaching S opt and then decreases. The critical temperature 
(proportional to A(T = 0) in the BCS approach) follows the 
same doping dependence. 

The superfluid density In the superconducting phase 
the thermodynamic properties at low temperature are es- 
sentially determined by the quasiparticles near the nodes, 
i.e. by the superconducting gap A, whereas the pseudo- 
gap amplitude A p plays no role. The superfluid density 
for instance can be evaluated according to the simple BCS 



0.035 
0.03 
0.025 
> 0.02 
j 0.015 
0.01 
0.005 




T=0 
T=8.0 meV 
T=8.5 meV 
T=9.5 meV 



5 10 15 20 25 30 35 40 45 
<p (degrees) 



Fig. 4. Temperature dependence of the LE for the same set 
of parameters of Fig. 3 (increasing temperature from top to 
bottom). At T > T c = 9 meV the superconducting gap closes, 
leaving a finite LEm ~ A p — 



formula: 



1 ^ f d 2 k }5% k / £„ k ,/?£„ k \ 



dh) dE^ J 



(8) 



where f(x) is the Fermi function. p s decreases linearly in T 
at low temperature, as expected in a d-wave superconduc- 
tor, with a slope determined only by the superconducting 
gap. For the assumed band dispersion (^|) it can be shown 
that 

T 16tln2 , , 

p a T)-p.(Q)~-- . (9) 

Zi 7T 

The slope of the superfluid density at low temperature, 
a = dp s /dT{T = 0) is estimated in Ref. 0] by means of 
direct ARPES measurements of the slope v/± of the super- 
conducting gap near the node. Because the observed va 
decreases in underdoped region (contrary to the gap at the 
M points), the doping dependence of a oc 1/va is found to 
increase by decreasing doping with respect to its optimum 
value. This behaviour is confirmed by our model: the slope 
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of dp s /dT given in Eq. (|^) is controlled by the supercon- 
ducting gap A, which follows the doping dependence of T c 
in the underdoped regime. However, as observed in Rcf. 
|0], the general trend of a(5), obtained by direct measure- 
ments of penetration depth, seems to be opposite, giving 
a decrease in the underdoped compounds of about 40% 
respect to its value at optimum doping. The correspon- 
dence between ARPES estimates and experimental data 
on a (S) requires the inclusion of doping dependent Lan- 
dau renormalization factors, as discussed in Ref. Q using 
the expression for p s with the inclusion of quasiparticles 
interaction (see also Ref. |n|). The presence of a dop- 
ing dependent Landau factor is plausible within our sce- 
nario if we attribute the origin of quasiparticle scattering 
to quasi-critical charge and spin fluctuations as in Eq. (|l|) . 
The same holds for a correspondence between the slope of 
the density of states iV(£) near £ = 0, where JV(£) rj ?y£ 
with 77 oc 1/va, and the values of rj extracted by specific 
heat measurements [Ira. 
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Fig. 5. p s (T) at different doping for A v = 0.16 eV, t = 0.25 
eV in weak coupling (V = 0.09 eV, left) and strong coupling 
(V = 0.9 eV, right). The temperature at which the numerical 
curves intersect the line 8T/ir is the Tkt, as explained in the 
text. Phase fluctuations are important between Tkt and Tbcs- 

related to the variation of the number of carriers with dop- 
ing: in a lattice model it does not increases monotonically, 
but reaches a maximum at an intermediate < 5 < 1 , and 
then decreases. Nevertheless, the doping for the maximum 
/9 S (0) does not coincides necessarily with 5 op t- 



The temperature dependence of the superfluid density 

As the temperature is increased, the system crosses n , . r ,-, r , 

J allows us to give an estimate ot the range of tempera- 

tures in which the phase fluctuations of the supercon- 



over to a regime of higher energy excitations, which sample 



regions where LE > Z\, and the slope is reduced. This can , , , , T , , , , 

' r ducting order parameter are relevant. Indeed, the criti- 

be seen in Fig. 5, where the temperature variation of pJT) r , , . ,. . , 

' y ' cal temperature tor phase coherence m a two-dimensional 



in weak and strong coupling is reported. The variation of 



system, Tkt, is related to the phase stiffness p s via the 



the slope with respect to its low temperature value is less T ^ , ,., r™ n , .. m / a i-u 

f f r Kosterhtz- 1 houless relation p s — 81 kt/tt- Assuming the 

pronounced in the strong coupling case, where almost the nf ,n , , , r m ,1 ■ 
^ o r ty 1 PCS temperature dependence for p s {l), the above re- 
same parameter controls the LEm and the slope of the n ,. , ,, Irn N 

r lation becomes a sell-consistency condition p s \±kt) = 

superconducting gap at the node. 0^1 j ^ i_ a ■ ^ 1 , 

00 oJkt/tt, and the phase fluctuations play an important 

Finally, we consider the doping dependence of the su- role in the range of temperatures between Tkt and Tbcs- 

perfluid density. The variation of the p s (0) with doping is As it is naturally, in our model this region is very narrow 



10 



L. Benfatto, S. Caprara, and C. Di Castro: Gap and pseudogap evolution 



in weak coupling, and becomes larger and larger as the 
strength of the pairing interaction increases. Indeed while 
p s (P) mainly depends on the doping 5, Tbcs rapidly in- 
creases with increasing the coupling, and becomes much 
larger than Tkt- In Ref. P2[ , due to the discontinuity 
of the density of states at the band edges, at low doping 
a regime T c ;g> Tkt is found. In such a case the transi- 
tion is expected to be Kosterlitz-Thouless. In our model 
the density of states at the Fermi level vanishes smoothly 
and continuosly as half-filling is approached. Therefore, 
in the weak-coupling limit, the critical temperature in- 
creases slowly with doping, and the Kosterlitz-Thouless 
regime is never found. Consistently with our assumptions 
within this model, in the regime of interest for cuprates, 
i.e. at weak coupling, the superconducting transition is 
essentially BCS-like, even in the underdoped region. The 
main part of the non classical behaviour is enforced by the 
presence of A p . 

Phase fluctuations effects have also been invoked to ex- 
plain the linear-T behaviour of p s (T) Within a model 
with pseudogap formation in the p-p channel Jll|] these ef- 
fects deserve a careful analysis. 

4 Conclusions and discussion 

In conclusion we want to summarize the previous results in 
closer connection with the phenomenology of the cuprates. 
The discussion of the previous section has given some pre- 
liminary indications on the set of parameters suitable to 
reproduce the general shape of the phase diagram ob- 
served in the cuprates. Nevertheless, to be consistent with 



the experiments, we have to introduce the temperature 
and doping dependence of the pseudogap parameter A p (T, S). 




Fig. 6. Left panel: doping dependence of T* , To and of the 
critical temperature T c . Right panel: doping dependence of the 
superconducting gap A(T = 0) and of the zero temperature 
LEm. The set of parameters used is t — 250 meV, V — 85 meV, 
c = 7, S = 0.2, and we assumed T (S) = 40 [1 - (5/S ) 4 } meV. 
The value of T*(5 = 0) = To (S — 0) has been extrapolated 
from Ref. §. 

The mechanism for pseudogap formation has not been 
the main issue in this paper, which focuses on the prop- 
erties of superconductivity arising in the pseudogap state. 
Nevertheless, according to the discussion of Sec. 1, we can 
imagine that A p {T, ^)7k schematizes the whole complica- 
tion of the strong scattering due to quasicritical charge 
fluctuations in the p-h channel, which arises in the stripe- 
QCP scenario ||. We assume that the doping dependence 
of A p follows the doping variation of the temperature To 
(To > T*) at which a first variation in the density of states 
occurs, as revealed by NMR and resistivity measurements 
M. We adopt the simple relation: 
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the optimum doping, as expected in the QCP scenario 
A (T S) — cTo(S)g (T/Tq(S)) (10) Finally, we interpolate between S — and Sq with the ex- 

pression To(<5) = 40[1 — (<5/<5o) 4 ] meV, which allows us to 

where g(l) = 0, g(0) = 1. g(x) interpolates smoothly be- 

reproduce the shape of the curve T*(S) observed in the 

tween these two limits ^ and c is a constant which we use 

underdoped regime. 

as a fitting parameter. Since the identification of T (S) 

via NMR and resistivity measurements leaves large un- The resulting phase diagram is shown in Fig. 6. Even 
certainty, we decide to refer to the ARPES experiments, within a simplified description of the doping and tem- 
which however only provide the temperature T* of LE perature dependence of A p , we recover the bifurcation 
closure. We proceed in the following way. In the heavily 01 T * and T c observed in the underdoped regime, while 
underdoped regime the main temperature dependence of th ey merge around optimum doping. T c follows the typi- 
the LE M (at fixed doping) is due to the decreasing of A p cal bell-shaped curve as a function of S. Our approach is 
with increasing the temperature, rather than to the tern- sti11 lacking of the temperature and doping dependence of 
perature variation of the chemical potential. Then T* and the superconducting coupling, which plays an important 
T coincide at low doping in our model, whereas the ex- role in determining the T C (S) variation around the QCP, 
perimental data seem to indicate that T stays larger then as it has been shown elsewhere [g|l0| . Even though the in- 



T*. This is an artifact of our approach, which is not de- traduction of this further complication would definitively 

voted to establish the connection between the pseudogap improve the agreement with the experimental data, with 

state and the Mott insulator and/or the antifcrromagnetic a faster deca y 01 T c at high doping, it would not change 

phase. This is a relevant open problem still under inves- the main features of the pseudogap state obtained here, 

tigation. Meanwhile we put T (S = 0) = T*(S = 0) and We reproduce the variation of the zero temperature LE M 

extrapolating to 5 = the T*(5) dependence reported in in all the phase diagram: according to the general dis- 

Ref. § for Bi2212 we obtain T (S = 0) = 40 meV. The cussion of the previous section, the LE M is determined 

constant c and the coupling V are adjusted to fix the val- b Y thc normal state pseudogap in the underdoped region, 

ues of the LE in the underdoped regime and the maximum and by the superconducting gap in the overdoped regime. 

T c . The temperature T is assumed to decrease with in- As a consequence, in the underdoped regime the LE M is 

creasing doping and to vanish at 5 = 0.2, slightly above uncorrected to the characteristic energy scale of low tem- 



2 Here we assume g{x) = (1 - x 4 /3)Vl - x 4 , which repro- perature quasiparticle excitations, which probe the value 

duces a mean-field-like behaviour near T = and To- However, of the superconducting order parameter around the nodal 

the specific form of g(x) does not play a crucial role in deter- points. Having now a doping and temperature dependent 

mining the general shape of the phase diagram. A p , the doping for the maximum T c does not coincide any- 
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more with the doping at which the normal state leading to Landau factors has been adressed. We observe here an- 
edge closes, as we had in the presence of a constant A p . other discrepancy: the slope of p s (T) near T c in our model 
Indeed, approaching optimum doping A p itself closes at a reduces by increasing the temperature contrary to the ex- 
lower temperature, and as a consequence the normal state perimental findings |Q. As already indicated in Section 
LEm closes faster than it would be under the only effect (3), within our model the superconducting fluctuations 
of the variation of the chemical potential with doping. are only important in the short region between Tkt aid 

Tbcs- The inclusion of these fluctuations, by reducing p s 



g = o.i2 near T c , would at least partially take care of this discrep- 
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Even though we did not reach yet the point to exploit 
fully the momentum and doping dependence of a quasi- 
singular effective interaction among quasiparticles arising 
nearby an instability in the p-p and p-h channels, alto- 
gether within the present simplified approach we produce 
a behaviour of T c (5), of T*(S) and of the LEm(0,S) in 
reasonable agreement with the experimental findings. 



Fig. 7. Doping dependence of the normalized superfluid den- 
sity p a (T) / p s (0) as a function of T/T c . We used the same set 
of parameters of Fig. 4. According to the phase diagram of 
Fig. 4, the optimum doping is at S opt w 0.16. Observe that in 

the overdoped regime S > 8 op t the superfluid density has the 

. j , , cijuu- We acknowledge A. Perali, C. Castellani, and M. Grilli 

standard d-wave mean field behaviour. b ' ' 
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